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Abstract. We use data on the local 3-dimensional galaxy 
distribution for studying the statistics of the detection 
rates of gravitational waves (GW) coming from supernova 
explosions. We consider both tensor and scalar gravita- 
tional waves which are possible in a wide range of relativis- 
tic and quantum gravity theories. We show that statistics 
of GW events as a function of sidereal time can be used for 
distinction between scalar and tensor gravitational waves 
because of the anisotropy of spatial galaxy distribution. 

For calculation of the expected amplitudes of GW sig- 
nals we use the values of the released GW energy, fre- 
quency and duration of GW pulse which are consistent 
with existing scenarios of SN core collapse. The ampli- 
tudes of the signals produced by Virgo and the Great At- 
tractor clusters of galaxies is expressed as a function of 
the sidereal time for resonant bar detectors operating now 
(IGEC) and for forthcoming laser interferometric detec- 
tors (VIRGO). Then, we calculate the expected number 
of GW events as a function of sidereal time produced by 
all the galaxies within 100 Mpc. 

In the case of axisymmetric rotational core collapse 
which radiates a GW energy of 10 -9 M Q c 2 , only the clos- 
est explosions can be detected. However, in the case of 
nonaxisymmetric supernova explosion, due to such phe- 
nomena as centrifugal hangup, bar and lump formation, 
the GW radiation could be as strong as that from a co- 
alescing neutron-star binary . For radiated GW energy 
higher than 10~ 6 Mqc 2 and sensitivity of detectors at the 
level h w 10~ 23 it is possible to detect Virgo cluster and 
Great Attractor, and hence to use the statistics of GW 
events for testing gravity theories. 

Key words: gravitation - relativity - waves - super- 
novae: general - galaxies: - clusters: general 
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1. Introduction 

In a few years the third generation of gravitational wave 
detectors will start searching for the most energetic events 
in the Universe caused by gravitational collapse and merg- 
ing of relativistic compact massive objects (see the re- 
view by Thorne 1997). This opens a new window onto the 
Universe and creates new connections between optical ex- 
tragalactic astronomy and gravitational wave astronomy. 
This will be the beginning of genuine observational study 
of the physics of the core collapse supernova explosions 
and testing relativistic and even quantum gravity theories 
(Damour 1999; Gasperini 1999). 

Expected sources of powerful gravitational wave (here- 
after, GW) events are connected with supernova explo- 
sions and merging of neutron stars and other relativis- 
tic compact massive objects in galaxies. Predicted GW 
signals essentially depend on the details of the last rela- 
tivistic stages of the gravitational collapse which is still 
poorly known (Thorne 1987;1997; Paczyncki 1999; Bur- 
rows 2000). Moreover studies of scalar-tensor gravity the- 
ories have shown that spherical gravitational collapse and 
binary systems generate scalar GWs which may be de- 
tected by existing GW detectors (Baryshev 1982; 1995; 
1997; Baryshevfe Sokolov 1984; Sokolov 1992; Damourfc 
Esposito-Fareze 1992; 1996; 1998; Shibata et al. 1994; Harada 
et al. 1997; Bianchi et al. 1998; Damour 1999; Brunetti et 
al. 1999; Gasperini 1999; Novakfe Ibanez 1999; Maggiore& 
Nicolis 2000 Nakao et al. 2000). 

The aim of this paper is to estimate the contribution 
of nearby galaxies and clusters of galaxies, within a radius 
of about 100/i^Q 1 Mpc around our Galaxy, to the detection 
of the possible GW events. This requires the knowledge 
of the actual 3-dimensional galaxy distribution (section 
2), the intrinsic rate of the most powerful events expected 
from different galaxy types (section 3), and the amplitude 
of the GW signal according to the prediction of existing 
scenarios of SN core collapses (section 4). In section 5 
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we calculate the probability of GW events as a function 
of sidereal time for currently operating bar detectors and 
forthcoming interferometric detectors. In this section we 
first study the amplitude expected for the Virgo cluster 
and the Great Attractor and then derive the density of 
probability of GW events as a function of the sidereal 
time for some detectors. A discussion of the results and 
the main conclusions are given in section 6. 

2. Distribution of galaxies within 100 Mpc 

Explosions of supernovae and mergings of binary massive 
compact objects are very rare in our Galaxy. Hence, only 
observations of many galaxies are expected to yield a rea- 
sonable detection rate. 

From the Lyon-Meudon extragalactic database LEDA 
we extracted a sample of 33557 nearby galaxies within 100 
Mpc. The 2D-distribution is shown on a Flamsteed equal 
area projection (Fig. |l|). Some prominent structures ap- 
pear. What are their distances? What is the actual space 
density of galaxies? 

2.1. The problem of the determination of distances 

The direct determination of distance requires difficult mea- 
surements and is thus available only for relatively small 
samples (5000 galaxies). The most efficient method con- 
sists of using the radial velocity together with a given 
Hubble constant. However, the determination of the Hub- 
ble constant is still controversial. Most of the difficulties 
come from the treatment of statistical biases. From the 
most recent absolute calibration and the best unbiased 
determination we adopt: 

Ho = 60 ± Skm.s^.Mpc- 1 (1) 

Nevertheless, we must consider that some results are 
still suggestive of a larger value (70 ± 5km. s^ 1 .Mpc" 1 ), 
and even H a = 50km. s _1 .Mpc -1 cannot be excluded. The 
smaller the value, the more difficult is GW detection due 
to larger distances to GW sources. 

2. 2. The problem of spatial distribution 

Analysis of the 3D-galaxy distribution from the correla- 
tion function method (Davis & Peebles, 1983; Davis, 1997) 
leads to the conclusion that the characteristic correlation 
length is Iq « 5Mpc and the maximum inhomogeneity 
scale Imax ~ IQMpc. However, a more general statistical 
method to study large-scale galaxy distribution has been 
recently developed (see the review by Sylos Labini, Mon- 
tuori, & Pietronero, 1998). It is applicable to any distri- 
bution of matter without the assumption of homogeneity 
(which is required in the correlation function analysis). 



The new analysis (the so-called conditional density func- 
tion approach) is actually taken from modern statistical 
physics where it works as a standard tool. 

Application of the conditional density function analy- 
sis to available redshift surveys of galaxies, such as CfA, 
SSRS, Perseus- Pisces, IRAS, LCRS etc, has revealed the 
fractal structure of the galaxy distribution up to the scales 
corresponding to the depth of these catalogs, i.e. about 
100 Mpc (see e.g. Pietronero et al.,1997; Sylos Labini et 
al.,1998). The fractal dimension of the spatial distribu- 
tion is close to T> w 2. From the KLUN galaxy survey 
(Teerikorpi et al., 1998) where distances to galaxies are 
obtained by the Tully-Fisher method, it was shown that 
the observed number-distance relationship corresponds to 
a fractal dimension T> « 2.2 and continues up to the depth 
of the KLUN catalog, i.e. 200 Mpc. 

The fractality implies that around any galaxy (includ- 
ing our own Galaxy) the density decreases as r v ~ z . This 
means that the number of galaxies does not increase as r 
but rather as r v ' . The direct consequence for the present 
analysis is that the detection rate of GW events will be 
lower than previously thought. 

From our sample we plotted the cumulative curves 
log AT vs. logr (N is the total number of observed galax- 
ies) within the radius r (Mpc) (Fig. ||). This is done for 
different absolute magnitudes. Intrinsically faint galax- 
ies (M = —17) start to be missed beyond logr = 1.3 
(«s 20Mpc), while galaxies brighter than M = —22 are 
observed up to the limit of our sample (r* lOOAfpc). The 
observed growth-curves correspond to V « 2.5 (dashed 
curve in Fig. ||). They are used to calculate the correc- 
tion allowing us to estimate the true number of galaxies 
in each direction, from the observed number. At a given 
distance, this correcting factor is simply deduced from the 
ratio of the observed and expected population (assuming 
T> rs 2.5). It is to be noted that up to logr = 1.3, even the 
faintest galaxies (M = —17) follow the linear curve. This 
means that the sample is complete up to this distance (20 
Mpc). 

3. Extragalactic sources of GW and intrinsic 
rates of events 

There are many sources of gravitational radiation in a 
galaxy. In fact, any accelerated motion generates GW. 
Among those usually discussed in the literature, galactic 
sources of GW are: supernova explosions, coalescing bi- 
nary systems, binary stars, rotating asymmetric pulsars, 
active galactic nuclei. We consider only the GW sources 
which arc expected to be sufficiently frequent and efficient 
to be detected in the near future. 

The most powerful sources of gravitational radiation 
are the core-collapse supernovae (types lb and II) , merging 
neutron stars (ns), and black holes or other relativistic 
compact massive objects (cmo) (Thorne 1987; 1997)) and 
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Fig. 1. Flamsteed equal area projection of our sample of 33557 galaxies located within 100 Mpc. 
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also supernovae of type la, which are probably due to the 
explosion of a CO white dwarf and might be sufficiently 
strong candidates. 

The relative supernova rates R g , s for galaxy type 'g' 
and SN type 's' are free parameters in our code. For further 
calculations we adopt these from van den Bergh & Tam- 
mann(1991), measured in SNU which equals one event per 
100 yr per 10 10 L^. In Table we give the adopted rates 
of SN events for different morphological types of galaxies. 
Note that we use the Hubble constant 60 km.s~ 1 Mpc~ 1 . 

The event rate R me rg for coalescing compact binaries 
composed of ns (or cmo) is still widely discussed and has 
a large uncertainty (e.g.Lipunov et al. 1997; Portegies- 
Zwart& McMillan 1999; Kalogera 1999). Here we adopt 
the values from Lipunov et al.(1995), however these events 
give a small contribution to the total statistics. 

The first detailed study of the gravitational wave sky 
produced by galaxies within 50 Mpc was done by Lipunov 
et al.(1995). They considered a wide class of GW sources 
in galaxies and used Tully's Nearby Galaxies Catalog com- 
prising 2367 galaxies. In this paper we use 33557 galaxies 
from the LEDA catalogue and consider both tensor and 
scalar GW from supernova explosions. 

4. Theoretical amplitudes expected for GW 
events from SN 



Table 1. Adopted intrinsic rate of GW events. We con- 
sider Supernovae explosions (SNIa, SNIb, SNII) and coa- 
lescences of White Dwarfs (WD), Neutron Stars (NS), and 
Black Holes or Compact Massive Objects (CMO). The 
rate may depend on the morphological type of the parent 
galaxy (Elliptical(E), Spirals(S) or Irregulars(Irr)). 



source Rate Rg lS of GW events 
per 100 years and 10 10 Lq 

SNIa 1.0 (E) 0.5 (S and Irr) 

SNIb 0.0 (E) 0.3 (S) 0.9 (Irr) 

SNII 0.0 (E) 1.4 (S) 4.2 (Irr) 
WD- WD 3 10~ 2 

NS-NS 3 10~ 3 

CMO-CMO 10" 3 



4-.1. The problem of supernova explosion 

Expected amplitudes and forms of GW signals from su- 
pernova explosions detected on the Earth by gravitational 
detectors essentially depend on the adopted scenario of 
core-collapsed explosion of massive stars and relativistic 
gravity theory. This is why the forthcoming GW astron- 
omy will give for the first time experimental limits on pos- 
sible theoretical models of gravitational collapse including 
the strong field regime and even quantum nature of the 
gravity force. 
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log(distance (Mpc)) 

Fig. 2. Cumulative curves drawn for a wide range of abso- 
lute magnitudes (from M = — 24 to M = — 17). The com- 
pleteness is severely affected for the less luminous galax- 
ies (M = —17) when the distance increases. The dashed 
curve shows the linear trend expected for a fractal dimen- 
sion 2.5. These curves are used to derive the true space 
density of galaxies. 

For the estimates of the energy, frequency and duration 
of supernova GW emission one needs a realistic theory of 
SN explosion which can explain the observed ejection of 
the massive envelope. Unfortunately, for the most inter- 
esting case of SNII explosion such a theory does not exist 
now. As was recently noted by Paczynski(1999), if there 
were no observations of SNII it would be impossible to 
predict them from first principles. 

Modern theories of the core collapse supernova are able 
to explain all stages of evolution of a massive star before 
and after the explosion. However, the theory of the explo- 
sion itself, which includes the relativistic stage of collapse 
where a relativistic gravity theory should be applied for 
the calculation of gravitational radiation, is still contro- 
versial and unable to explain the mechanism by which the 
accretion shock is revitalized into a supernova explosion 
(see the discussion by Paczynski 1999; Burrows 2000). 

Moreover, recent observations of the polarization of 
core collapse supernovae (Wang et al. 1999) and the rela- 
tivistic jet in SN1987A (Nisensonfe Papaliolios 1999) give 
strong evidence in favor of a jet-induced explosion mecha- 
nism for massive supernovae (Khokhlov et al. 1999; Mac- 
Fadyen & Woosley 1999; Wheeler et al. 2000). Further evi- 
dence for highly asymmetric SN explosions comes from re- 
cent observations of afterglows and host galaxies of gamma- 
ray bursts (Paczynsky 1999). This means that new "non- 
standard" scenarios of SN explosions (and hence GW ra- 
diation) may appear in the future and it may become im- 
portant to study different possibilities for the expected 
GW signal (hence a wide range of GW parameters). Here 
we adopt scenarios existing in the literature to estimate 



the GW signal and we postpone the discussion of non- 
standard possibilities to section 6. 

4-2. Standard pulses of gravitational radiation 

The main aim of the present paper was not to provide real- 
istic supernova explosion models, but to study the statis- 
tics of GW signals which may be expected in standard 
SN explosion scenarios. Hence we do not enter into the 
detailed calculations of the precise forms of GW signals 
within different gravity theories, but will simply use gen- 
eral energy arguments. 

In our calculations we use the standard pulse of gravi- 
tational radiation introduced by Amaldife Pizzella(1979), 
which is a GW burst of sinusoidal wave with amplitude 
ho, frequency vo and duration r g . For the case of tensor 
GW, the amplitude ho of the signal on the Earth due to 
the GW burst that occurs at a distance r with total energy 
E gw is (see e.g. Pizzella 1989): 

h Q = 1.4 x l0- 2 °(^)(^)(-^)i/2(i£)i/2 (2) 
r u IMqC 2 T g 

For the case of scalar GW the amplitude ho is given 
by the same relation with a pre-factor depending on the 
considered theory, e.g. for tensor field gravity theory it is 
w 2 x 10~ 20 (Baryshev 1995). 

Hence, each type of GW source at a fixed distance r is 
characterized by three main observable parameters E gw , 
vo and r g . In the next subsection we choose acceptable 
values for these parameters. 

4-3. Adopted values for energy, frequency and duration 
of tensor and scalar GW pulses 

Let us first consider the parameters for standard ten- 
sor GW pulses in General Relativity. There is no unique 
widely accepted model for tensor GW radiation produced 
by SN core collapse and in the literature two main scenar- 
ios are usually discussed: axisymmetric and nonaxisym- 
metric ones, reviewed by Thorne (1997). 

Within the theory of axisymmetric rotational core col- 
lapse, Zwerger & Muller (1997) found that the energy 
spectrum covers a frequency 50 Hz < v Q < 3 kHz but 
most of the power is emitted between 500 Hz and 1 kHz. 
Duration of the pulses lies between 0.5 — 5 ms. According 
to numerical calculations by Stark & Piran (1985), con- 
sidered in Ferrari et al. (1999) as a basis for prediction of 
GW background from SN, the GW energy spectrum has 
two maxima around 5 kHz and 9 kHz. The duration of 
GW pulses also is of the order a few ms. In accordance 
with these calculations, for our statistical approach, we 
adopt the characteristic frequency v Q = 1 kHz and the 
duration of the pulse r 9 = 1 ms. 

For the total GW energy radiated by SN core col- 
lapse there is a very large range of predicted values in 
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the literature. According to Zwerger & Muller (1997) the 
energy radiated in the form of GWs lies in the range 
6. lfJ- u M c 2 < E gw < 8. 10~ 8 Af Q c 2 , which is consis- 
tent with Bonazzola & Marck's (1993) results for the de- 
formation parameter s < 0.1. However fully relativistic 
numerical simulations by Stark & Piran, also adopted by 
Ferrari et al. (1999), give E gw < l(T 3 M Q c 2 . Moreover, 
if the collapsing core rotates so rapidly that it becomes 
nonaxisymmetric and may be transformed into a bar-like 
configuration, which also might break up into several frag- 
ments, then the GW radiation could be almost as strong 
as that from a coalescing neutron star binary. Several spe- 
cific scenarios for such nonaxisymmetric SN core collapses 
have been proposed (see review by Thorne,1997). Accord- 
ing to Lai & Shapiro (1995) the energy radiated in GW 
during the nonaxisymmetric stage of the gravitational col- 
lapse can be as large as 4 x 10~ 3 M Q c 2 . Bonncll & Pringle 
(1995) considered gravitational radiation from SN core 
collapse and fragmentation, which could produce the GW 
energy 10 _2 M Q c 2 . For our statistical study of GW events 
we choose, as a basis, the value E gw = 10 _6 M Q c 2 . 

Within classical general relativity there is no scalar 
GW and for instance spherically symmetric collapse does 
not generate gravitational radiation. However other rel- 
ativistic and quantum gravity theories predict both ten- 
sor and scalar GW. Calculations of amplitudes, frequen- 
cies and forms of the scalar gravitational radiation in the 
case of spherically symmetric SN core collapse has been 
made by Shibata et al.(1994), Harada et al.(1997), No- 
vak & Ibanez (1999). The released scalar GW energy is of 
order ui~ 1 Mqc 2 , where u> > 10 3 is the parameter of the 
Brans-Dicke theory. The characteristic frequency and du- 
ration are similar to the tensor GW in general relativity. 
For the comparison of the statistics of tensor and scalar 
events we adopt the same energy, frequency and duration 
for both scalar and tensor GWs. 

In Fig.|| we plot the expected amplitudes of GWs cal- 
culated according to Eq.2 for a wide range of the parame- 
ter E gw , which covered the expected values of GW energy 
for tensor and scalar GW. Three pairs of lines (a,b,c) cor- 
respond to the following combinations of the main GW 
pulse parameters: 

a) E gw = 10~ 3 M©c 2 ; v = 10 3 Hz; r g = lms; 

b) E gw = l(r 6 M Q c 2 ; v = 10 3 Hz; T g = lms; 

c) E gw = lCT 9 M Q c 2 ; v = 10 3 Hz; r g = lms. 

For further calculations we adopt the values of GW 
amplitudes which correspond to the case b. Obviously, our 
calculations may be rescaled by using any other combina- 
tions of main GW parameters according to Eq.2. 

In Figure ^ we draw two levels of sensitivity (horizontal 
dotted lines) which can be expected today (ho = 1CT 21 
and 10~ 22 ). One sees that it would be possible to detect 
SN explosions at the distance of the Virgo cluster and of 
the Great Attractor only with an optimistic energy release 
(case a). 
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Fig. 3. Theoretical GW amplitude versus distance. The 
predicted amplitude at a distance r is given for tensor 
(solid lines) and scalar (dashed curves) waves according to 
Eq.2 for cases a,b,c (see text). For instance, case b corre- 
sponds to the GW energy of 10 -6 Af Q c 2 , a pulse duration 
T g = 1CP 3 s and a frequency i/q — 1 kHz. Two sensitiv- 
ity levels are represented (dotted lines) corresponding to 
ho = 1CT 21 and 1(T 22 . The distances of LMC, M81, Virgo 
and Great Attractor are labelled on the x-axis. 

4-4- Longitudinal and transversal scalar GW 

An important difference between tensor and scalar GWs 
is that tensor waves (spin 2) are transversal while scalar 
waves (spin 0) can be transversal and/or longitudinal. 
There is no scalar GW in general relativity. 

In the frame of the Jordan-Fierz-Brans-Dicke theory, 
as for any metric tensor-scalar theories, the scalar wave 
is transversal but isotropic in the plane transversal to the 
propagation direction (Damour & Esposito-Farese (1992). 
The transversal spin wave may be presented as the met- 
ric perturbation in the following form (Bianchi et al. 1998; 
Brunetti et al. 1999; Maggiore & Nicolis 2000): 

9ik = ho cos(cji — kz)diag(0, 1, 1, 0) (3) 

In the frame of the field approach to gravity formulated 
in 60's by Feynman (see e.g. Feynman et al. 1995; Thirring 
1961; Straumann 2000), there are also scalar gravitational 
waves which are longitudinal and generated by the trace of 
the energy-momentum tensor (Baryshev 1982,1995; Sokolov 
1992). In such a scalar wave, a test particle moves along 
the direction of the wave propagation. The most straight- 
forward way to demonstrate this is to use exact relativistic 
equations of motion of test particles in tensor and scalar 
potentials derived by Kalman (1961). It is an important 
fact that the physical interaction of a GW with a detec- 
tor may be completely analyzed in terms of the weak field 
approximation, i.e. with the usual Lagrangian formalism 
of relativistic field theory in Minkowski space (see Ap- 
pendix) . 
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The scalar plane monochromatic GW in the system of 
coordinates with the z-axis directed along the wave prop- 
agation may be presented in the form: 

hik = h(t, z)r) ik (4) 

where h(t, z) is a 4-scalar field of Minkowski space, or 

Kk = h cos(u>t — kz)diag{\ 1 — 1, — 1, —1) (5) 

where h$ <C c 2 is the amplitude of the wave, u> = 2iri/Q, k = 
uj/c — 2ir/X, rjik — diag(l, — 1, — 1, — 1) is the Minkowski 
metric tensor. As shown in the Appendix the test par- 
ticles in such a gravitational potential are moving in the 
direction of the wave propagation (z-axis). The analysis of 
the motion of test particles is sufficient for the description 
of an interaction of the scalar GW with interferometric 
and bar detectors because the scalar GW does not inter- 
act with electromagnetic field. Indeed, the interaction La- 
grangian is A int = fkkT^ = due to the tracelessness of 
the energy-momentum tensor of the electromagnetic field. 
Hence without any ambiguity, the scalar GW is longitu- 
dinal in the field gravity theory. 

This means that both types of scalar GWs, longitudi- 
nal and transversal, are theoretically possible and phys- 
ically different. The physical difference between longitu- 
dinal and transversal scalar GW may be experimentally 
established by comparing the direction of the maximum 
sensitivity of a bar detector with the direction of the axis 
of the bar. Indeed, for a bar detector the maximum sensi- 
tivity to the transversal GW is in the direction orthogonal 
to the bar axis, while for the longitudinal scalar GW the 
maximum sensitivity direction is along the axis of the bar. 

For an interferometric detector (which has two arms) 
the maximum sensitivity to the tensor GW will be in the 
direction orthogonal to the plane containing both arms. 
For the longitudinal scalar waves there are two directions 
of maximum sensitivity which coincide with the directions 
of each arm. It is interesting that for the transversal scalar 
GW the directions of maximum sensitivity also coincide 
with the two arms of the interferometer (Maggiore,Nicolis 
2000; Nakao et al.2000). This is a special case of a two 
arm interferometer, while in the case of a bar detector 
the detector patterns are different for longitudinal and 
transversal scalar GWs. 

The difference in the response of GW detectors to ar- 
riving tensor and scalar GW pulses allows one to test the 
nature of the detected waves. Below we compare the statis- 
tics of expected GW events for transversal and longitudi- 
nal GWs. 

5. Probability of GW events as a function of 
sidereal time 

5.1. The expression for the number of GW events 



be either detected or not, depending on the amplitude 
of the signal and on the sensitivity of the detector. The 
number of detected events produced at a time t for a given 
type of wave (tensor or scalar noted by index w) and a 
given detector (interferometer or bar detector noted by 
index d) will be the sum of the contribution of each galaxy 
to the considered type of source (SNIa, SNIb, SNII etc... 
noted by index s). 

— (i, s, to, d) — L g x R g>s x A (6) 
g 

where L g is the luminosity (in L Q ) of the g-th galaxy; R is 
the rate of GW events per year and per Lq for the consid- 
ered type of source and the considered type of wave ; A is 
the detectability calculated from the observed amplitude 
h bs and the limiting amplitude hu m of the detector. 
Let us detail the three terms used in Eq. ||. 

• The luminosity L g is derived from the absolute blue 
magnitude M g of the g-th galaxy. A correction for in- 
completeness is done using Fig|2| and assuming that 
missed galaxies are distributed like the observed ones. 
This does not account for structures completely hid- 
den by the disk of our Milky Way 0. The luminosity 
is: 

L g = C(M g , r9 )iQ-O-4(M 9 -M ) (?) 

(r g is the distance of the g-th galaxy and M Q = 5.2 is 
the absolute magnitude of the Sun in the photometric 
B-band) From Fig.^j it is visible that the correcting 
term C is almost equal to one below 20Mpc (logr = 
1.3), because all curves are linear even for the faintest 
galaxies. This means that up to the Virgo cluster (ss 
18Mpc) all galaxies are included. 

• The rate R 9iS depends on the morphological type of 
the considered galaxy (index g) and on the type of the 
GW source (index s) according to Table |l|. As an ex- 
ample, we made calculations for SNII + SNIb. Hence, 
the expected counts will be underestimated, if all other 
conditions are satisfied. 

• The detectability A is defined as: 

^ J 1 if h Q bs ^ hn 7n , ,gx 

\ if hobs < h Um . 

We will consider a sensitivity of hu m = 10 -23 . It is not 
yet accessible by present and scheduled detectors, but 
one can hope an improvement will be possible from 
space missions. The observed amplitude is defined as: 

h obs = h .\G(a, S,w,d)\ (9) 

h is given by Eq.||. It depends on the considered type 
of source (via the energy, frequency and duration of 
the pulse) and on the considered type of gravitational 
wave. G(a, 8, w, d) is the geometrical factor. It charac- 
terizes the relative orientation of the GW and of the 



Let us consider a sample of galaxies. Each galaxy may 
produce GW events with a certain rate. These events will 



1 For example, the Great Attractor which is presented in 
section 5.2 
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Fig. 4. The main geometrical definitions. Z is the zenith 
of the site. P is the northern pole of the equatorial coordi- 
nate system. 7 defines the sidereal time with respect to the 
southern meridian. The source S is located by its equa- 
torial coordinates a = jOH for the right ascension and 
5 = HOS for the declination. The distance to the zenith 
is ( = ZOS. The reference direction for the detector is the 
direction OX. It is either the x-arm for an interferometric 
detector or the direction of the bar for a bar detector. The 
azimuth of this direction is $ D . It is defined in the direct 
sense from the north to the west. 

detector at the time of the observation. It depends on 
the type and position of the detector (characterized 
by its latitude, sidereal time and azimuth of its refer- 
ence axis), on the direction of the GW (characterized 
by the equatorial coordinates a, S of the g-th galaxy), 
and on the type of wave (including polarization for 
tensor waves). Note that the longitude of the site is 
not needed explicitly because it is included in the def- 
inition of the sidereal time. The sidereal time will be 
given in hours (from to 24). 

The calculation of the geometrical factor G is explained 
in the next subsection for interferometric and bar detec- 
tors. 



5. 2. Geometry of the system 

The geometrical configuration is defined by the wave direc- 
tion (OS), the local horizon (perpendicular to OZ , where 
Z is the zenith) and the equatorial plane (perpendicular to 
OP, where P is the equatorial northern pole). The main 
features of the system are shown in Fig.^. 

The expression of the geometric factor for tensor waves 
has been calculated by Schutz & Tinto (1987) and Thorn 
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(1987). For scalar waves it has been calculated by Bary- 
shev (1997). 

The relevant angles are £, $, O and VP. They depend 
on the different cases. Let us detail : 

• For interferometric detectors: the zenith distance £ = 
ZOS (angle between the direction of the zenith and the 
arrival-direction of the wave), the azimuth $ = XOK 
of the direction of the wave (direction OK) calculated 
with respect to the reference axis of the_detector. 

• For bar detectors: only the angle 9 = XOS is relevant. 

• For tensor GW the polarization angle 'F As in Schutz 
& Tinto (1987) it may be measured in the GW plane 
from the line of nodes (intersection of the GW plane 
and of the horizontal plane) . 

Note that the direction of the polarization depends 
on the geometry of the emitting source. For the calcula- 
tion with the actual galaxy distribution \& will be simply 
chosen at random. With these notations we obtain the 
following results: C an d $ are calculated by solving the 
spherical triangle ZPS. One obtains: 

cos £ = sin L sin S + cos L cos 5 cos(i — a) (10) 

L is the latitude of the site where the detector operates. 
a and (5 are the equatorial coordinates of the source, and 
t is the sidereal time at the site. £ is defined in the range 

(0,7T). 

cos L sin (5 — sinLcos<5cos(t — a) ,,,, 
cos($ + $) = ^ '- (11) 

y/1 — COS z Q 

cos<5sin(i — a) . . 

sin($ + $) = 1 > (12) 

y 1 — cos^ (, 

$ is the azimuth of the reference axis of the system 
counted in the direct sense from the direction of the north. 
In this paper the reference axis is OX, as in Thorn 1987. 
It is either the direction of the X-arm for an interferomet- 
ric detector or the direction of the bar for a bar detector. 
$ is defined over the range (0, 2n). 

The angle is calculated by solving the triangle XZS: 

cos 6 = sin £ cos 4> (13) 

O is defined over the range (0,7r). 

The relevant angles being calculated according to the 
previous relations, the geometrical factor in the four con- 
sidered cases has the following form: 

1. For tensor waves and interferometric detector the 
geometrical factor is: 

G = 0.5(1 + cos 2 C) cos 2$ sin 2* - cos ( sin 2$ cos 2* (14) 

2. For longitudinal scalar waves and interferometric 
detector the geometrical factor is: 

G = sinC(| cos$| - I sin<f>|) (15) 

3. For tensor waves and bar detector the geometrical 
factor is: 

G = sin 2 6 cos 2* (16) 
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4. For longitudinal scalar waves and bar detector the 
geometrical factor is (Appendix, Eq. A22): 



We considered the following detectors: 



G = cos e 



(17) 



The interferometric GW detectors (such as TAMA, 
GEO600, VIRGO, LIGO) have the frequency range of 
about 10- 1000-Hz and the sensitivity h 1(T 20 - 1CT 22 . 
Sensitivity is a measure of the detectable amplitude and 
it is proportional to Al/l, the relative length variation of 
the arms of the interferometer. 

Presently working bar detectors, such as IGEC RBO 
(International Gravitational Event Collaboration of Reso- 
nant Bar Observatory) which includes five cryogenic reso- 
nant bar detectors (ALLEGRO, AURIGA, EXPLORER, 
NAUTILUS, NIOBE), have a typical bandwidth of the 
order of 1 Hz around each one of the two resonances 
(close to 1000 Hz). The achieved sensitivity is now Fo ~ 
2 - 6 x lO" 21 ^- 1 (Prodi et al. 2000). In future one ex- 
pects the sensitivity to be at the level of 10~ 22 Hz~ 1 . 

We have calculated the predicted amplitude for dif- 
ferent existing detectors and specific galaxy clusters. The 
closest concentration of galaxies is the Virgo cluster. The 
name of the VIRGO detector comes from the name of the 
cluster itself because it may be the main source of first 
detectable GW events. The Virgo cluster is not very far 
from our Local Group (20 Mpc). It induces a velocity of 
about 170 km.s^ 1 on the Local Group. On the other hand, 
it has been claimed (Dressier et al., 1987) that a hidden 
large concentration of galaxies (hereafter, the Great At- 
tractor) induces a velocity of about 500km. s^ 1 on our 
Local Group (i.e. about three times more). The distance 
of such a concentration has been estimated to be three 
times the distance to Virgo. This means that the number 
of galaxies could be about 27 times larger than the num- 
ber of Virgo galaxies 0. If the sensitivity of GW detectors 
is improved, the Great Attractor may become the major 
source of GW events. This justifies the interest we place 
in this region. The position of this putative Great Attrac- 
tor would be roughly at galactic coordinates I = 310 deg, 
b = Odeg. This is well supported by the apparent 2D- 
distribution of galaxies which shows that this region may 
constitute a link between two visible structures on both 
sides of the Milky way (Paturel et al., 1987 and Fig.|l|) and 
by the discovery of a large number of galaxies around this 
region (Kraan-Korteweg, 2000). 

Then, we considered two clusters as dominant sources. 
The adopted equatorial coordinates and distances are the 
following: 



Cluster 



a(1950) de/to(1950) 



Virgo 

Great- Attractor 



12/i28m 
15/iOOm 



-12 deg 40' 
-60 deg 00' 



20Mpc 
GOMpc 



2 The velocity, and thus the acceleration during the same 
time, is three times larger. The distance is also three times 
larger, resulting in a mass nearly 3/3 -2 = 27 times larger 



Detector 


Latitude L 


Azimuth <I> 


Type 


VIRGO 


44.7 deg 


— 15 deg 


intcrf. 


AURIGA 


45.4 deg 


+45 deg 


bar 


NAUTILUS 


41.8 deg 


Odeg 


bar 


NAUTILUS 


41.8 deg 


90 deg 


bar 



In the following subsections we calculate the ampli- 
tudes in different conditions, for different detectors and 
sources. For the study of the polarization effect on ten- 
sor waves (Fig.|| to ||) we used the cluster Virgo (as a 
point source) and detectors VIRGO and AURIGA. For 
the comparison of the distribution of the amplitudes of 
Virgo and of the Great- Attractor (Fig. [?] to [l(]) we used 
Virgo and the Great attractor as point sources and the 
detectors VIRGO, AURIGA and NAUTILUS. For the cal- 
culation of the density of probability of GW events along 
the sideral time (Fig. [IT] and [12]) we used the sample of 
individual galaxies and the detectors VIRGO, AURIGA 
and NAUTILUS. 



5.3. Effect of the polarization 

Because the polarization angle cannot be predicted, it is 
important to show the effect of the polarization for tensor 
waves. Considering only the Virgo cluster and interfero- 
metric (VIRGO) and bar detector (AURIGA), we calcu- 
lated the amplitude according to Eq. || using h a for case b 
(Fig.||). The geometrical factor is given by Eq. [[J and ||, 
for an interferometric and bar detector respectively. We 
used 36 polarization angles over the range (0, 2ir) The 
results are shown in Fig. |^a-b , respectively. 

For interferometric detectors, the different curves are 
shifted along the x-axis (sidereal time) depending on the 
polarization. On the contrary, for bar detectors, the curves 
are shifted along the y-axis (amplitude), but the minima 
and maxima always appear at the same x values (same 
sidereal times). Let us explain why it is not correct to 
calculate the mean over the different polarizations. If an 
event is produced with the favorable polarization, it will 
be detected if the amplitude is larger than the limiting 
amplitude. On the other hand, if the polarization is unfa- 
vorable, the observed amplitude will be reduced and the 
considered event may fall below the limiting amplitude. 
The global effect of the uncertainty on the polarization is 
simply a reduction of the number of counted events, but 
not a reduction of the amplitude. Finally, the distribution 
of the amplitudes along the sidereal time will be simply 
given by the envelope of the curves obtained for the differ- 
ent polarizations. It must be noted that it is a handicap for 



3 In fact, the period is 7r/2 for interferometric detectors and 
7r/4 for bar detectors 
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sideral time (VIRGO) 



5 10 15 

sideral time (VIRGO) 




sideral time (AURIGA) 

Fig. 5. a) Amplitude as a function of sidereal time for ten- 
sor GW emitted by Virgo and seen with the VIRGO inter- 
ferometric detector for different polarizations of the GW. 
Each curve represents 36 fixed polarization directions, b) 
Amplitude as a function of sidereal time for tensor GW 
emitted by Virgo and seen with the AURIGA bar detector 
for different polarizations. 

interferometric detectors because the contrast is smoothed 
and the total number of expected events is reduced. We 
repeated the same calculations with polarizations taken 
at random (Fig. ||a-b). The same effect is clearly visible. 
The curves are the envelopes of the previous ones. Some 
events are seen with smaller amplitudes due to unfavor- 
able polarization. 



5.4- Comparison of the amplitudes of GW from Virgo 
and the Great Attractor 

We have calculated the GW amplitudes as a function of 
sidereal time for the four detectors listed in the previous 
table. For each of them we give two figures, for tensor 



sideral time (AURIGA) 

Fig. 6. a) Amplitude for the VIRGO detector as in Fig.||- 
a but with random directions of polarization. The shape 
is unchanged but some events have a reduced amplitude. 
This will lead to a reduction of the GW events detected 
at a given sensitivity level, b) Amplitude for the AURIGA 
detector as in Fig.|5|b but with random directions of po- 
larization. The shape is unchanged but, on average, the 
amplitude is reduced. 



and scalar GW, respectively. In each figure we present the 
amplitudes expected for sources in the Virgo cluster (solid 
curve) and in the Great Attractor region (dashed curve). 

The results are given in four figures (Fig. @ to §3): The 
caption of each figure contains detailed comments. Here, 
we will simply highlight prominent features. 

With a sensitivity hu m = 10 -23 only GW events from 
the Virgo cluster (solid lines) will be detectable. The Great 
Attractor (dashed curves) will be detectable only with a 
sensitivity of hu m = 10~ 235 . Nevertheless, there is one 
case (VIRGO detector - Fig.^b) where scalar waves could 
be seen from the Great Attractor with a sensitivity of 
Because we expect about 25 times more 



10 



-23 
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events from the Great Attractor, this may result in a peak 
in the rate of events around sidereal time t = 17h. 

Another important features in these diagrams of am- 
plitudes is that tensor and scalar waves give peaks which 
generally have opposite phases. In other words, if we ex- 
pect a maximum of events for scalar waves, there should 
be a minimum for tensor waves. This is clearly visible by 
comparing Fig.|| a and b. This is also important because it 
may be used to disentangle the contributions of these two 
kinds of waves. The same characteristic is present when 
we compare the expected counts for the NAUTILUS bar 
detector with two perpendicular orientations (<I> = deg 
and $ = 90 deg). The figure suggests that NAUTILUS 
could benefit from an orientation complementary to the 
one used, e.g., with the AURIGA orientation. 

5.5. Number of GW events expected for the actual galaxy 
distribution within 100 Mpc 

Using the catalog of galaxies described in section 2, we 
simulated the count of GW events using Eq. 0. The calcu- 
lation is made with an energy of E gw — 10~°Mqc 2 , a fre- 
quency of 1 kHz and a duration of 1 ms. This corresponds 
to the case b of Fig|| The adopted limiting amplitude hn m 
is 1CU 23 . It should allow us to reach the Virgo cluster. The 
calculation is done for VIRGO, AURIGA, NAUTILUS 
(azimuth=0 deg) and NAUTILUS (azimuth=90 deg). We 
plotted simultaneously the counts for tensor (solid curves) 
and scalar waves (dashed curves). The results are shown 
in Fig.|ll] and [l2| for the four considered detectors. When 
comparing these Figures with the expected amplitudes 
calculated for the detectors VIRGO and AURIGA, it is 
seen that all predicted maxima are at the expected po- 
sitions according to the distribution of amplitudes given 
by the Virgo cluster alone (solid curves in Fig. |7^ and b 
and ||a and b). This confirms that, within a distance of 
20 Mpc, the Virgo cluster dominates. This was not obvi- 
ous because the influence of other galaxies was not easy to 
predict. The expected number of GW events reaches a few 
tens per year at the most favorable sidereal time but with a 
yet unreachable sensitivity for the considered GW energy. 
For instance, from Fig.|Tl|-a we see that several maxima 
are expected for scalar waves (dotted curve) in agreement 
with their positions predicted from Fig.^-b for the Virgo 
cluster alone. Similarly, the maximum at t ks 13/i for ten- 
sor waves (solid curve) is the one predicted for Virgo alone 
as seen from Fig.^-a. This confirms that the Virgo clus- 
ter will be the dominant source of GW events when the 
sensitivity will be hu m = 10 -23 . 

From Fig.0-b (AURIGA bar detector) one can see 
clearly that tensor waves and scalar waves have opposite 
phases as far as the sidereal time is concerned. The same 
effect is also visible from Fig.[l2] a and b with the NAU- 
TILUS bar detector. Further, we see that changing the 
orientation of the bar by 90 deg also produces a change 
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Fig. 7. a) Amplitude as a function of sidereal time for 
tensor GW emitted by Virgo (solid curve) and the Great 
Attractor (dashed curve) as seen with the interferometric 
VIRGO detector. With a sensitivity of 10~ 23 only Virgo 
will give large enough amplitude to make events detectable 
between t = 9h and t = 15h and between 23ft- and 4h. 
The Great Attractor could be detected with a sensitiv- 
ity of 10 -23 - 5 (except between t = llh and t = I5h). b) 
Ibidem but for scalar waves. Virgo (solid curve) could be 
detectable at different peaks along the sidereal time. The 
Great Attractor could be barely detected with a sensitiv- 
ity of 10~ 23 at t = 17 h 

of sidereal time phase; The maximum in Fig.|l2]-b for, say, 
tensor waves (solid curve), corresponds to a minimum in 
Fig-H-a. 

6. Discussion and conclusions 

From the detection of GW events with bar and interfer- 
ometric detectors, one can see that, due to the different 
geometrical factors and to the anisotropy of the distribu- 
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Fig. 8. a) Amplitude as a function of sidereal time for 
tensor GW emitted by Virgo (solid curve) and the Great 
Attractor (dotted line) as seen with the bar detector AU- 
RIGA. With a sensitivity of 1CT 23 Virgo will be detectable 
between t = Oh and t = 18h. b) Ibidem for scalar waves. 
Virgo could be detectable at two main positions t = 9h 
and t = 2lh. The Great Attractor could be barely detected 
with a sensitivity of 10 -23 at t = 12h. 
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Fig. 9. a) Amplitude as a function of sidereal time for 
tensor GW emitted by Virgo (solid curve) and the Great 
Attractor (dashed curve) as seen with the bar detec- 
tor NAUTILUS oriented in the direction of the north 
(azimuth=0 deg). This figure is the same as Fig.^-a but 
with a shift in sidereal time, b) Ibidem with NAUTILUS 
oriented in the east- west direction (azimuth=90 deg). 
Virgo could be detectable at two main positions t = Ih 
and t = 12h. Figures a and b have opposite phases. 



tion of the GW sources, it is in principle possible to make 
the distinction between transversal tensor GW, transver- 
sal scalar GW and longitudinal scalar GW. 

To demonstrate this difference we calculated the ex- 
pected amplitudes and the number of events as a func- 
tion of sidereal time for real positions of existing GW de- 
tectors and for real 3-dimensional distribution of galax- 
ies within 100 Mpc. If one adopts the value for the en- 
ergy of GW pulse E gw = 10~ 6 M Q c 2 (lines b in Fig.|), 
then, GW events produced at the distance of the Virgo 
cluster can be detected only with a sensitivity, yet un- 
reachable, of hu m = 10 -22 5 . However if the GW energy 
is about E gw = 1O _3 M0C 2 (case a in Fig. 3) as predicted 



by nonaxisymmetric scenarios of SN core collapse, then 
Virgo cluster and Great Attractor would be visible with 

hlirn = 10~ 22 . 

If the GW energy emission has value E gw = 1O _9 M0C 2 
or less, as predicted by an axisymmetric scenario of SN ex- 
plosion, then Virgo would simply be not detectable with 
present and forthcoming detectors. Thus, one cannot ex- 
pect a high detection rate, but only serendipitous detec- 
tions from very nearby SN's. 

We would like to emphasize a point which seems im- 
portant to us. Today there are several scenarios of GW ra- 
diation from SN core collapse even in the frame of the Gen- 
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Fig. 10. a) Amplitude as a function of sidereal time for 
scalar GW emitted by Virgo (solid curve) and the Great 
Attractor (dashed curve) as seen with the bar detec- 
tor NAUTILUS oriented in the direction of the north 
(azimuth=0 deg). This figure is the same as Fig.^-a but 
with a shift in sidereal time, b) Ibidem with NAUTILUS 
oriented in the east-west direction (azimuth=90 deg). 
Virgo could be detectable at two main positions t = 7h 
and t = 18h. Again, Figures a and b have opposite phases. 



eral Relativity, which predicts only tensor gravitational ra- 
diation (see review by Thorne,1997). Within different sce- 
narios, the SN core collapse may lead to a large range of ra- 
diated GW energy from lO _11 M c 2 up to lO~ 2 M c 2 and 
to very different forms of GW bursts, i.e. to the different 
spectral energy distributions and durations of GW pulses 
from msec up to sec and even minutes (Lai & Shapiro, 
1995). It is important to recall that other relativistic and 
quantum gravity theories (such as string theory, the Jordan- 
Fierz-Brans-Dicke theory and tensor field theory) predict 
scalar gravitational radiation which is generated also in 
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Fig. 11. a) Density of probability of GW events (num- 
ber per year) for the VIRGO interferometric detector as a 
function of the sidereal time of the site. The calculation is 
made from the actual sample of galaxies described in sec- 
tion 2. Tensor waves are shown with a solid curve, while 
scalar waves are given as dashed curves. These predicted 
counts are obtained from Eq.^ using case b of Fig|| for 
the calculation of the observed amplitude. The sensitivity 
is assumed to be hu m = 1CV 23 . Because one retrieves the 
expected distribution found for Virgo alone (solid curves 
of Fig.0), one concludes that the Virgo cluster will be the 
dominant source of GW events with such a sensitivity, b) 
Ibidem for the AURIGA bar detector. The comparison 
with predicted amplitude for Virgo cluster alone also con- 
firms that it will be the dominant source of GW events at 
the considered sensitivity. It is clearly visible that tensor 
and scalar waves have opposite phases. 
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the case of a spherical gravitational collapse and which 
may carry a large GW energy. 

We would like to emphasize the importance of con- 
sidering a wide range of GW burst parameters for SN 
core collapses by discussing the SN1987A and SN1993J 
events. Analysis of data recorded by Geograv for SN1987A 
(Amaldi et al. 1987) and by Explorer- Allegro for SN1993J 
(Mauceli et al. 1997), showed that there are GW candidate 
events, which the authors themselves do not consider as 
real signals because the GW energy calculated for a stan- 
dard pulse with a duration r g — Ims (and hence band- 
width 1 kHz) gives E gw « lO 3 M c 2 for both supernovae. 
However, in the case of pulse duration of about 1 sec (and 
hence bandwidth 1 Hz) the GW energy needed for pro- 
ducing the same GW amplitude is about 1 Mqc 2 . In this 
case the observed GW amplitudes correspond to the lines 
"a" in Fig. 3 and fit well the decrease expected from the 
relative distances of the two host galaxies (LMC, M81) 0. 

This means that GW data analysis should be done 
for the interval of possible signal durations from ms to 
sec timescales. For the GW pulses longer than 1 sec the 
frequency bandwidth is less than 1 Hz and the sensitivity 
of bar detectors may be essentially improved if one uses 
as a signal the difference between two signals coming from 
two resonances of a bar detector. 

Let us conclude with the most secure results presented 
in this paper: 

— We wrote a computer code to calculate the amplitude 
of transversal or longitudinal gravitational waves from 
any real source which could be detected with bar or 
interferometer detectors. This code allows us to fix all 
experimental conditions: latitude and sideral time of 
the detector, type and orientation of the detector, right 
ascension and declination of the emitting source, char- 
acteristics of the gravitational wave. 

— This code has been used to calculate how the ampli- 
tude changes with sideral time of the site for existing 
detectors. In particular, using a point source we study 
the effect of the polarization of the wave in the case of 
tensor waves. 

— Then, we compared the distribution of amplitudes along 
the sideral time for bar and interferometric detectors 
and for transversal and longitudinal gravitational waves 
for two dominant point sources. We found an interest- 
ing result for bar detectors (NAUTILUS) which can 
have different orientations. This may be used for the 
selection of real GW signals from the noise. The result 
is that it would be theoretically possible to discrim- 
inate both kinds of waves because, for a considered 

4 The high value of energy release in a gravitational collapse 
could also be supported by the existence of gamma-ray bursts 
if they are related to SN explosions because they may generate 
in the electromagnetic channel an energy of about 1 Mqc 2 with 
timescale a few seconds (Kulkarni et al. 1999; Paczynski 1999). 



point source, the maxima of the distributions do not 
appear at the same sideral time. 
— Finally, we applied the code to calculate the expected 
count of events generated by the actual distribution 
of galaxies using the gravitational energy release pre- 
dicted by existing scenarios of Supernova core collapse. 
The result is that, with the sensitivity of GW de- 
tectors h — 1CP 22 and for the released GW energy 
E gw = 1CP 3 , predicted by nonaxisymmetric scenarios 
of SN core collapse supernova, it is possible to detect 
GW events from the distance of the Virgo cluster and 
Great Attractor and to use the statistics of the events 
as a test of gravity theories. 

Appendix A: The longitudinal character of the 
scalar wave in the field gravity theory 

Feynman's field gravity theory (Feynman et al.,1995) is 
based on the Lagrangian formalism of the relativistic quan- 
tum field theory and presents a non-geometrical descrip- 
tion of gravitational interaction. According to Feynman, 
the gravity force between two masses is caused by the 
exchange of gravitons which arc mediators of the gravita- 
tional interaction and actually represent the quantum of 
the relativistic tensor field ip lk in Minkowski space rf k . 

It is important that the problem of the physical in- 
teraction of a gravitational wave with an detector may be 
completely analyzed in terms of the weak field approxima- 
tion where classical relativistic field theory is applicable 
(see Landau & Lifshitz, 1971). 

First, let us consider the general problem of the mo- 
tion of a relativistic test particle having rest-mass uiq, 
4- velocity u l , and 3- velocity v in the gravitational field de- 
scribed by the symmetric tensor potential ip lk in the flat 
Minkowski space-time. The Cartesian coordinates always 
exist and the metric tensor is rf k — diag(l, — 1, — 1, — 1) 
(we utilize notations of the text-book Landau & Lifshitz, 
1971 ). 

To derive the equation of motion for test particles in 
the frame of the field gravity theory we start from the 
stationary action principle in the form 

5S = <5(~ J (A (p) + A M) )dO) = (Al) 

where the variation of the action is made with respect to 
the particle trajectories 5x l for fixed gravitational poten- 
tial ip lk (t,r), and d£l is the element of the 4- volume. 
The free particle Lagrangian is 

A ( P ) = -V*T$) (A2) 
and the interaction Lagrangian is 

A (mt) = -^ipikTfa (A3) 

where the energy-momentum tensor of the test point par- 
ticle is 

Tfc = m Q c 2 S(r - r p ){l - W (A4) 
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The result of variation gives the equations of motion in 
the form (Kalman 1961; Baryshev 1986): 



dzt k 

A\— = -m Bl lU k u l 



(A5) 



where p k — niocu k is 4-momentum of the test particle, 



Al = (1 - ^i> nl u n u l )ri 



2 2 

—ip kn u n u l + —i> l k 



1 / ■< 



(A6) 



(A7) 



In Eq.(A5) the rest mass of the test particle may be can- 
celed, hence in the field gravity theory m.;„ = m g = too 
without the initial equi vale nce postulate (Thirring 1961). 
In agreement with Eq.(A5) all classical relativistic post- 
Newtonian gravity effects have the same values as in gen- 
eral relativity (see e.g. Baryshev 1995). 

The scalar gravitational wave in the field gravity the- 
ory is generated by the trace of the energy momentum 
tensor T — 7]ikT tk of the sources of the gravitational po- 
tential (Baryshev 1982;1995; Sokolov 1992) and may be 
expressed in the form: 

iP ik (t,x) = ^(t,x)r, ik (A8) 

where ip = rjikip lk is the trace of the tensor potential ip lk 
and hence is a 4-scalar field in Mink owski space. 

Substituting Eq.(A8) in Eq.(A5) and taking into ac- 
count the weak field approximation we get the following 
equation of motion of the test particle in scalar gravita- 
tional potential: 

% -a* <"> 

Spatial components of this equation (i — a) give the 
expression for the gravity force acting on the test particle: 

dp ™- -. 2 

dt 



m n . , v" , , v . dtp v , « , , . 

[(l-_)VV> + ~(^ + ~-W)](A10) 



V 1 - 

The corresponding 3-acceleration of the test particle is 

dv _ 1 
~dt ~ ~ 



v dtp . 



(All) 



The time component (i = 0) of the Eq.(AS) gives the work 



produced by the scalar gravitational wave: 
dEk mo t ^ ! i v 2 dtp , 

~df 



-.{v ■ Vip 



4J1 



c 2 cdt 



(A12) 



where is the kinetic energy of the particle. For the case 
of slow motion of the test particle (v/c <C 1) we get for 
the gravity force: 
dp mo, 

~dt ~ 4" 

for 3-acceleration of the test particle: 

| = - Jv* (A14, 



-Vtp 



(A13) 



and for the work of the gravity force: 
dE k m 

— = — — (v-Vtp) (A15) 

The scalar plane monochromatic GW in the system of 
coordinates with the z-axis directed along the wave prop- 
agation may be written in the form: 

ipik = ji>(t, z) rj ik (A16) 
ipik = ho cos(ujt — kz) diag(l, —1, —1, —1) (A17) 

where ho <C c 2 is the amplitude of the wave, o> = 2irv, 
k = uj/c = 2n/X. 

For the scalar gravitational potential Eq.( A17| ) the 3- 
acceleration Eq.(Al4) of the test particle is 

^ = ^ = (A18) 
dt dt y ' 

dv z 1 dip 

~dT ~ Tcdt 

According to Eq.(Al9) under the action of the scalar grav- 
itational wave a test mass undergo small longitudinal os- 
cillations in the direction of wave propagation (z-axis): 

Xjio 
2tt c 2 

The relative distance between two particles Az(t) = z\ (i) — 
z 2(t) = l(t) gives the relative oscillations Al(t) = l(t) — 
lo = AIq cos(ujt + a), which may be detected by gravita- 
tional detector. 

The very important property of the scalar gravita- 
tional potential is that it does not interact with the elec- 
tromagnetic field. Indeed, the interaction Lagrangian for 
the potential Eq.([Alj|) is 
1 



(A19) 



z(t) 



-sinlut — kz$) + zq 



(A20) 



A 



(int) 



(em) 



= o 



(A21) 



because the trace of the energy momentum tensor of the 
electromagnetic field T( em ) = 0. This means that the above 
analysis is sufficient for the study of the response of the bar 
or interferometric detectors where the distances between 
test particles are measured by means of the electromag- 
netic field. 

For a detector with the length lo <C A the amplitude 
of the relative displacement A^o / lo between two test par- 
ticles, which is measured by the detector is 

Ajo 



ho cost 



(A22) 



where ho is the amplitude of the scalar wave Eq.(Al7), 
and is the angle between the direction of the wave prop- 
agation and the line connecting the two masses. 

For a bar detector the integral cross section for the 
scalar gravitational wave having Av g > Av a is given by 
(Baryshev, 1997): 



^(scalar) 



= ^^MV s 2 cosH 

7T cr 



(A23) 



where M is the mass of the cylinder, V s = loqL/tt is the 
speed of sound in the cylinder, loq is resonance angular 



Baryshev Yu.V., Paturel G.: Rates of gravitational waves 



15 



frequency of the bar, L is the length of the cylinder and 
6 is the angle between the bar axis and the direction of 
the wave propagation. The factor cosO in Eqs.( A22 A23 ) 
shows the longitudinal character of the scalar gravitational 
wave. 
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sideral time for NAUTILUS (azimuth=90deg) 

Fig. 12. a) Same as the previous figure for the NAUTILUS 
pointing towards the north (<E> = deg). The results 
and the conclusions are very similar to those obtained 
from Fig.[Tl[- b. Nevertheless we note that this orienta- 
tion is less favorable for scalar waves in comparison with 
the AURIGA bar detector, b) Ibidem for the NAUTILUS 
pointing towards the west ($ Q = 90 deg). We note the im- 
provement of detected GW events when comparing with 
the previous orientation. 



